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•^he  problem  of  a balanced,  planar  or  eucisymmetric , supersonic  jet 
impingi-ig  normally  on  a flat  surface  has  been  considered  based  on 
an  inviscid  theory.  The  object  of  the  study  was  to  provide  a 
rational  model  for  calculating  shock-interference  heating  as 
produced  by  a type  IV  shock- interaction  pattern.  The  unwanted 
singularity  at  a low  supersonic  Mach  number  peculiar  to  scheme  I 
of  the  one-strip  foinnulation  of  the  method  of  integral  relations, 
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s Observed  by  South  and  by  Guimner  and  Hunt/  was  successfully 
removed  by  the  application  of  the  scheme  III  of  the  one-strip 
formulation  of  the  method  of  integral  relations.  The  resulting 
simultaneous  nonlinear  algebraic  equations  were  easily  solved 
iteratively  by  the  Newton- Raphson  method.  Sensitivity  of  the 
solution  on  various  approximating  functions  employed  was 
extensively  investigated.  Unlike  the  findings  reported  by  Gummer 
and  Hunt,  solutions  that  satisfy  all  well-posed  boundary 
conditions  can  be  obtained  by  the  one-strip  formulation.  Results 
indicate  that,  for  the  planar  case,  a rational  engineering 
solution  for  the  stagnation-point  velocity  gradient  (and  hence  the 
peak  heat- transfer  rate)  has  been  obtained.  For  the  axisymmetric 
case,  however,  solutions  appear  to  be  not  quite  converging.  A 
two-strip  formulation  based  on  the  method  of  integral  relations 
is  also  included. 
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NORMAL  IMPINGEMENT  OF  A SUPERSONIC  JET  ON  A PLANE  - A BASIC 
STUDY  OF  SHOCK- INTERFERENCE  HEATING 


This  report  presents  a theoretical  method  to  predict  the  severity 
of  shock- interference  heating  caused  by  the  impingement  of  a shock 
wave  on  a blunt  fin.  The  problem  of  a supersonic  jet  (resulting 
from  the  interaction  of  the  incident  shock  with  the  fin  bow  shock) 
impinging  on  the  fin  surface  was  studied  based  on  the  one-strip 
formulation  of  the  method  of  integral  relations.  A rational 
engineering  rolution  for  the  stagnation-point  velocity  gradient 
(and  hence  the  peak  heat-transfer  rate)  has  been  obtained  for  the 
planar  case.  The  present  jet-impingement  model  could  be  coupled 
with  the  shock- interference  model  of  Edney  to  predict  type  IV  shock- 
interaction  effects. 

The  present  study  was  sponsored  by  the  Naval  Air  Systems  Command , 
AIR-320C,  under  Air  Task  No.  A320-320C/WR023-02-003. 
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SYMBOLS 


speed  of  sound 


specific  entropy  function,  p/p' 

equal  to  zero  (planar  case)  or  one  (axisynunetric  case) 

Mach  number 
static  pressure 

2 2 1/2 

total  speed,  (u  + v ) 
coordinate  axis  along  the  plate  surface 
velocity  component  in  the  r-direction 
velocity  component  in  the  y-direction 
free-stream  velocity  of  the  jet 

coordinate  axxs  perpendicular  to  the  plate  surface 

constant,  iy  - l)/2y 

ratio  of  (constant)  specific  heats 

the  angle  the  upper  boundary  of  the  wall  jet  makes  with 
respect  to  the  negative  y-direction  (see  Fig.  1) 

detachment  distance  of  the  shock  wave  or  of  the  wall-jet 
boundary  (see  Fig.  1) 

location  of  the  sonic  point  at  the  wall 

the  angle  the  flow  behind  the  shock  wave  makes  with  respect 
to  the  negative  y-direction 

density 

the  angle  the  shock  wave  makes  with  respect  to  the  negative 
y-direction  (see  Fig.  1) 

method  that  employs  the  equation  of  global  mass 
conservation,  Eq.  (41) 

method  that  employs  the  equation  of  modified  continuity, 

Eq.  (7) 

method  that  employs  piecewise  smooth  approximating  functions 
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SYMBOLS  (Cont'd) 

SP  method  that  imposes  the  condition  of  satisfying  Eq.  (56) 

at  r = 0 

Subscripts 

j at  the  upper  boundary  of  the  wall  jet 

s at  the  shock  wave 

w at  the  plate  surface 

n at  the  surface  sonic  point 

0 at  r = 0 

1 at  the  line  of  the  jet  edge,  r = 1 

2 at  r = 1/2 

0“  at  free  stream 


INTRODUCTION 

As  an  extraneous  shock  wave  impinges  on  a blunt  body  in  a 
hypersonic  flow,  greatly  increased  aerodynamic  heating  and  pressure 
over  a very  small  region  near  the  impingement  point  have  been 
observed  (Refs.  (1)  to  (5)).  The  incident  shock  wave  may  be 
generated  either  by  boundary-layer  separation  (Refs.  (3)  to  (5))  or 


(1)  Edney,  B. , "Anomalous  Heat  Transfer  and  Pressure  Distributions  on 
^lunt  Bodies  at  Hypersonic  Speeds  in  the  Presence  of  an  Impinging 
Jhc-k,"  xFA  Report  115,  The  Aeronautical  Research  Institute  of 
Sweden,  Stockholm,  1968 

(2)  Hainfc , F.  D.  and  Keyes,  J.  W. , "Shock  Interference  Heating  in 
Hypersonic  Flows,"  AlAA  Journal,  Vol.  10,  1972,  pp.  1441-1447 

(3)  Hiers,  R.  S.  and  Loubsky,  V7.  J. , "Effects  of  Shock-Wave  Impinge- 
ment on  the  Heat  Transfer  on  a Cylindrical  Leading  Edge,"  NASA 
TN  D-3859,  Ames  Research  Center,  Moffett  Field,  Calif.,  1967 

(4)  Kaufman,  L.  G. , III,  Korkegi,  R.  H.  and  Morton,  L.  C,,  "Shock 
Impingement  Caused  by  Boundary  Layer  Separation  Ahead  of  Blunt 
Fins,"  ARL  TR  72-0118,  Aerospace  Research  Laboratories,  WPAFB, 
Ohio,  1972 


(5)  Giilerlain,  J.  D. , Jr.,  "Experimental  Investigation  of  a Fin- 
Cone  Interference  Flov/  Field  at  Mach  5,"  NSWC/WOL/TR  75-63 
Naval  Surface  Weapons  Center,  White  Oak  Lab.,  Silver  Spiing,  Md., 
1976 
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by  an  extraneous  surface  (Refs.  (1)  to  (3)).  Six  different  types  of 
shock-interaction  patterns  have  been  classified  by  Edney  based  on  an 
extensive  experimental  study  (Ref.  (1)).  Among  them,  the  type  IV 
interference  pattern . produces  the  most  severe  shock-in' .•'rference 
heating  and  pressure.  This  interference  results  in  a supersonic  jet 
embedded  in  the  subsonic  flow  field.  In  fact,  peak  interference 
heating  rates  up  to  17  times  the  interference-free  stagnation-point 
value  and  peak  pz'essures  up  to  eight  times  the  free-stream  pitot 
pressure  level  have  been  measured  by  Hains  and  Key^.s  (Ref.  (2)). 

Despite  its  significance,  past  analyses  (Refs.  (1)  to  (3)  and 

(6))  on  the  type  IV  interference  were  inadequate  and  generally 
empirical  in  nature.  Recently,  a time-dependent  finite-difference 
method  was  used  by  Tannehill,  Holst  and  Rakich  (Ref.  (7))  to  solve 
the  Navier-Stokes  equations  for  the  two-dimensional  shock- 
impingement  problem.  Although,  in  principle,  their  computer  program 
can  be  used  to  compute  all  six  types  of  shock  interactions,  only 
type  III  interference  results  have  been  published  so  far.  However, 
the  elaborate  computations  involved  and  the  extensive  computer  time 
required  by  their  method  make  it  highly  desirable  to  have  some 
relatively  simple,  yet  reasonably  accurate,  approximate  method. 

Such  an  approach  has  in  fact  been  pursued  by  Edney  (Ref.  (1))  and  by 
Keyes  and  Hains  (Ref.  (6)).  However,  their  empirical  treatments  of 
the  jet- impingement  process  suggest  the  need  for  a more  rational 
study.  This  is  the  subject  of  the  present  paper. 

The  impingement  of  a balanced  supersonic  jet  on  a flat  surface 
was  studied  both  theoretically  and  experimentally  for  an  axisymmetric 
jet  at  normal  impingement  by  Gummer  and  Hvmt  (Ref.  (8)),  and 
theoretically  for  a plane  jet  at  an  arbitrary  angle  with  the  surface 
by  Bukovshin  and  Shestova  (Ref.  (9)).  Both  groups  have  used  the 
scheme  I of  the  method  of  integral  relations  in  its  crudest  form 
(one  strip)  (Ref.  (10)).  However,  in  both  studies  the  centered 

(6)  Keyes,  J.  W.  and  Hains,  F.  D.,  "Analytical  and  Experimental 
Studies  on  Shock  Interference  Heating  in  Hypersonic  Flows," 

NASA  TN  D-7139,  Langley  Research  Center,  Hampton,  Va. ^ 1973 

(7)  Tannehill,  J.  C.,  Holst,  T.  L.  and  Rakich,  J.  V.,  "Numerical 
Computation  of  Two-Dimens ionaJ  Viscous  Blunt  Body  Flows  with 
an  Impinging  Shock,"  AIAA  Paper  75-134,  AIAA  13th  Aerospace 
Sciences  Meeting,  20-22  Jan  1975 

(8)  Gummer,  J.  H.  and  Hunt,  B.  L. , "The  Impingement  of  a Uniform, 
Axisymmetric,  Supersonic  Jet  on  a Perpendicular  Plat  Plate," 

The  Aeronautical  Quarterly,  Vol.  XXII,  Part  4,  1971,  pp.  403-420 

(9)  Bukovshin,  V.  G.  and  Shestova,  N.  P.,  "Incidence  of  Plane 
Supersonic  Jet  on  a Plane  at  an  Arbitrary  Angle,"  Fluid  Dynamics, 
Vol.  2,  No,  4,  1967,  pp.  97-100 

(10)  Belotserkovskii,  0.  M. , ed.,  "Supersonic  Gas  Flow  Around  Blunt 
Bodies,"  NASA  Technical  Translation  TTF-453,  Cjne  1967 
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expansion  to  3mb?-nt  pressure  of  the  jet-edge  streamline  behind 
the  shock  v/ave  vdz  not  properly  considered,  and  instead  an  empirical 
condition  of  sonic  velocity  at  the  jet  edge  behind  the  shock  wave 
was  imposed.  Turthermore,  at  low  supersonic  Mach  numbers,  both 
Sou  til  (Ref.  (11);  and  Gummer  and  Hunt  (Ref.  (8))  have  pointed  out 
the  singular  behavior  of  the  governing  equation  of  the  scheme  I of 
the  method  of  integral  relations.  This  singularity,  which  has  no 
counterpart  in  an  exact  solution,  will  cause  the  computation  in  the 
shock  layer  to  break  down.  This  is  of  special  importance  to  us 
since,  according  to  Edney  (Ref.  (1)),  low  supersonic  Mach  numbers 
are  in  the  range  of  particular  interest  to  the  shock-interference 
problem. 

The  singular:  ' can  be  shown  to  be  easily  removed  if  the 

governing  differential  equations  are  integrated  .nee  again  along  the 
body-surface  direction.  This  constitutes  the  scheme  III  of  the 
method  of  integral  relations  (Ref.  (10)).  This  approach  was 
utilized  in  the  present  study  to  generate  solutions  to  the  one-strip 
approximation  equations  of  the  jet-impinge-aent  problem.  As  we  shall 
shov.'  later,  in  contrast  to  th>  findings  reported  by  Gummer  and 
Hun.'  (Ref  (8)),  the  one-strip  approximation  does  yield  solutions 
the  . ..atii  fy  >11  well -posed  boundary  conditions.  A two-strip 
fo'  ation  of  the  problem  h-s  also  been  completed,  but  solutions 

hav»  'jt  yet  t.  .an  carried  out.  For  the  sake  of  completeness,  this 

is  i.  .'ads'"  i?'  v. 'u"  ^ppendi■t. 

PROBLEM  FORMULATION 


GOVERNING  EQUATIONS 

Cons  icier  the  flow  geometry  schematically  shown  in  Figure  1. 

The  origin  of  the  coordinate  system  is  placed  at  the  stagnation 
point  of  the  flat  surface.  The  problem  is  considered  to  be  steady 
and  two-dimenu lonal  or  axisymmetric,  with  r and  y axes  along  and 
perpendicular  to  the  plate  surface,  respect-vely , and  the  free- 
stream  jet  flow  is  in  the  negative  y-direction.  For  simplicity,  the 
gas  is  assumed  to  be  inviscid  and  obeys  the  perfect  gas  law;  its 
conditions  are  characterized  by  the  pressure,  p,  density,  p, 
temper?ture,  T,  and  velocity  components,  u and  v,  in  the  r and  y 
directions,  respectively.  Ahead  of  the  shock  wave,  the  jet  is 
assumed  to  be  uniform  with  constant  static  pressure  equal  to  the 
ambient  value.  These  assumptions  are  of  the  usual  kind  that  are 
generally  made  by  other  investigators.  Heat-transfer  rates  can  be 
calculated  using  the  well-known  boundary-layer  results  once  the 
pressure  distribution  along  the  plate  surface  is  determ’’  ed  from  the 
inviscid  approach. 


(11)  South,  J.  C.,  Jr.,  "Calculation  of  Axisymmetric  Supersonic 

Plow  Past  Blunt  Bodies  wxth  Sonic  Corners,  Including  a Program 
Description  and  Listing,"  NASA  TN  D-4563,  Langley  Research 
Center,  Hampton,  Va.,  1968 
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Under  these  conditions,  the  governing  conservation  equations 


|j:(r^pu)  + ~(r^pv)  = 0 


Ijr(r^puv)  + |p-[r^(Bp  + pv^)]  = 0 


~[r^(6p  + pu^)]  + |-(r^puv)  = jSp 


p =:  p(l  - q^) 


where 


2 = (Y  - 1? 

2y  ■ 

2 2 2 
q = u + V 

j = 0 or  1 for  two-dimensional  or  axisymmetric  jets,  respectively j 
and  Y is  the  ratio  of  (constant)  specific  heats.  The  variables 
are  all  nondimens ional.  Thermodynamic  variables  are  non- 
dime-nsionalized  by  the  corresponding  stagnation  values  in  the  free- 
stream  jet,  velocities  by  the  maximum  adiabatic  velocity  and 
distance  by  the  jet  radius.  Obviously,  the  magnitude  of  the  non- 
dimensional  free-stream  jet  velocity  is  related  to  the  free-stream 
jet  Mach  number  by  , 


= 

00 


(Y  - dm: 

2 + (Y  - DM? 


There  is  also  a geometric  relation 


in  the  shock  layer,  and 


= - cot  o 


- - cot  6 


in  the  wall-jet  layer,  where  e is  the  detachment  distance  of  the 
shock  wa'v's  or  of  the  wall  jet,  a and  6 are  the  angles  the  shock 
wave  and  the  upper  boundary  of  the  wall  jet  make  with  respect  to 
the  free-streaun  jet  flow  direction,  respectively  (see  Pig.  1) . 

The  method  of  integral  relationF  requires  that  the  governing 
partial  differential  equations  be  cast  into  divergence  form,  such 
as  Equations  (1)  to  (3) . However,  combinations  of  these  equations 
can  also  be  reprasented  in  divergence  form.  For  example,  one  may 
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combine  the  relation  of  constant  entropy  along  streamlines,  the 
energy  equation  (4),  and  the  continuity  equation  (1)  to  yield  a 
modified  continuity  equation 


- r 

3r  L 


r^ud 


o 1/(Y-1) 


1/{Y-1) 


0 


(7) 


which  was  the  original,  widely  employed  formulation  of 
Belotserkovskii  (Ref.  (12)).  For  a sphere  in  supersonic  flow, 

Xerikos  and  Anderson  (Ref.  (13))  found  that  the  one-strip  formulation 
based  on  the  modified  continuity  equation  yielded  results  which 
agree  with  experiments  better  than  that  based  on  the  original 
continuity  equation.  The  difference  is  expected  to  disappear  when 
the  number  of  strips  increases.  In  the  present  one-strip  formu- 
lations, however.  Equation  (7)  will  be  used  instead  of  Equation  (1). 


An  additional  simplification  arises  when  only  one  strip  is 
used  in  the  formulation,  namely,  the  strip  boundaries  are  either 
the  shock  wave  or  streamlines.  Along  the  plate  surface,  the 
constant  entropy  relationship  can  be  used  to  relate  pressure  to  the 
surface  velocity.  This  algebraic  relation  can  thus  be  employed  to 
replace  the  radial  momentum  equation  (3),  as  we  shall  see  in  the 
next  section. 


The  flow  field  can  be  divided  into  two  regions,  a shock-layer 
region  (0  - r ^ 1)  and  a wall-jet  region  (1  £ r ^ ri),  where  r = n 
is  the  location  of  the  sonic  point  at  the  wall 

u^(n)  = a^(n) 

and  it  is  unknown,  a priori.  The  two  regions  are  related  by  the 
requirements  that,  at  r = 1,  e,  E and  are  continuous  and  a and  6 
are  governed  by  the  Prandtl-Meyer  expansion  relation,  where  E is 
the  specific  entropy  function 


E = p/p^ 


and  <ii  is  the  stream  function.  If  0 is  the  angle  the  flew  behind  the 
shoc)c  wave  makes  with  respect  to  the  negative  y-direction,  then  the 
oblique  snock  relations  give 


cot  02^  = 


(y  * 

- 1) 


tan 


(8) 


(12)  Belotserkovskii,  0.  M. , "Flow  With  a Detached  Shock  Wave  About 
a Symmetrical  Profile,"  Journal  of  Applied  Mathematics  and 
Mechanics,  Vol.  22,  1958,  pp.  279-296 


(13)  Xerikos,  J.  and  Anderson,  W.  A.,  "An  Experimental  Investigation 
of  the  Shock  Layer  Surrounding  a Sphere  in  Supersonic  Flow," 
AIAA  Journal,  Vol.  3,  1965,  pp.  451-457 
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where  the  subscript  1 denotes  quantities  evaluated  at  r = 1.  Now, 
5^  is  related  to  by 


- tan  - 1)  I 


where 


2 . ^<3j 

“j  (y  - i)Pj  (Y  - 1)  (1  - q2) 

and  , 

2 2q2 

“si  ==  ^ T~ 

(Y  - 1)  (1  - q‘i) 

The  subscripts  j and  s denote,  respectively,  quamtities  evaluated 
at  the  upper  boundary  of  the  wall  jet  and  right  behind  the  shock 
wave.  Obviously, 

Pj  = ?»  = ] (3 


pj  = 


E.  = E , 
D si 


g.  - (1  - p./Pj)^'^  (13 

The  specific  entropy  function  evaluated  right  behind  the  shock  at 
r = 1,  depends  only  on  y and  Hence,  from  Equations  (8) 

to  (13) , we  obtain 

» fun(M„,Y,aj^) 


Since  the  upper  boundary  of  the  wall-jet  layer,  y = e (r)  for  r i 1, 
is  a streamline,  we  have 
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de  _ 4-  ^ 

^ . cot  6 


(14) 


first  two  quadrants  only,  we  may  combine  Equations  (4) 
and  (14)  to  yield 


Uj  = q^  sin  6 


Vj  = -qj  cc  6 


(15) 

(16) 


The  signs  are  determined  from  the  fact  that  u . - 0 for  r ^ 1 The 
boundary  conditions  are:  ^ ~ 


A.  At  the  wall,  y = 0 


= 0 

^w  " ®s0 


(17) 

(18) 


where  is  the  specific  entropy  function  evaluated  right  behind 
the  shock  at  r •=  0. 

B.  At  the  centerline,  r = 0 


= 0 

(19) 

II 

n 

CD 

O 

(20) 

- ir/2 

(21) 

1 the  shock  wave,  y = e (r) , r - 1,  the  Remkine-Hugoniot 

relations  for  the  gas  apply: 


® L(Y  + 


(M^sin^ 


dm! 


CT  - l)j 


(22) 


“ v.f 


2(M£sin^g  - 1) 


■[ 


1 + ~ dm; 


(Y  + DM, 
-1/(Y-D 


(y  t DM^sin^g 
,2  + (y  - Dm  sin^g 


] 


(23) 


(24) 


10 
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2YM£3in^g  ~ (y  - 1) 


ry  + 1) 


- l)tCsin‘‘g 

DM^sin^a 
'00  • 


p = E p’ 

D.  At  the  jet  boundary,  y = e(r),  r ^ 1,  Equations  (8)  to  (16) 
apply . 

METHOD  OF  INTEGRAL  RELATIONS  - SCHEME  I 

A.  SHOCK-LAYER  REGION.  Integrating  the  axial  momenttim 
equation  from  0 to  e,  and  utilizing  the  identity  .that 


e(r) 


e(r) 


/ |7(r='puv)dy  = Ip  / r^puvdy  - gfr^PgU^v 


we  obtain 


3T } • * • 

r^puvdy  - ifr^PjUsV^  + r^(S(Pg  - P„)  + - p„v^)  = 0 (?T, 


In  the  first  approximation,  the  integrand  is  assumed  to  be  linear  in 
y so  that  Equation  (27)  is  approximated  by 

~ 2r^{6(Pg  - p^)  + PgVg{Vg  + UgCot  a) } = 0 (28) 

Equations  (5)  and  (17)  have  been  used  in  the  above  equation. 
Similarly,  Equation  (7)  can  be  integrated  over  the  thickness  of  the 
shock  layer  to  yield 


I r^cj^u 


„ l/(y-l)  5 1/(Y-1) 

^a-q2,  +>s,a-u2) 


. 5 1/(Y-1) 

+ 2r-'(l  - q„)  [v„  + u„cot  a]  = 0 

s s s 

From  Equations  (4)  and  (18)  and  the  definition  of  the  specific 
entropy  function,  we  obtain  the  algebraic  relation  that 


2 y T 1/  (y~^  ) 

1 

^SO 


Since,  for  fixed  values  of  M^^  and  y,  the  quantities  evaluated  at  the 
shock  depend  only  on  o (as  can  be  seen  from  the  Rankine-Hugoniot 


» ■ 
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relations).  Equations  (5),  (28)  and  (29)  are  the  governing  equations 
for  the  variables  e,  o and  u . This  constitutes  the  scheme  I of  the 
method  of  integral  relations;  Initial  conditions  are  Equations  (19) 
and  (21)  . It  is  well  known  in  related  blunt-body  problems  that  idle 
missing  third  initial  condition  is  supplied  by  the  regularity 
condition  at  the  surface  sonic  point  (Refs.  (10)  to  (i2)).  For  the 
jet- impingement  problem,  this  requires  the  consideration  of  the  wall 
jet  since  the  surface  sonic  point  lies  outside  the  shock  layer 
(Ref.  (8)).  Before  we  proceed  any  further,  it  is  important  to  point 
out  a singular  feature  of  the  scheme  I formulation.  The  singularity 
occurs  as 


do 


0 


in  Equation  (28)  and  ^ becomes  unbounded.  This  has  no  counterpart 

in  an  exact  solution.  As  was  remarked  by  South  (Ref.  (11))  and  by 
Gummer  and  Hunt  (Ref.  (8)),  the  singularity  occurs  in  the  shock 
layer  for  - 2.  In  fact,  Gummer  and  Hunt  found  no  solution  that 

‘3(P«u  V ) 

will  satisfy  the  wall-jet  relations.  Since  will  appear 

in  any  method  that  approximates  the  integral  in  Equation  (27)  by  an 
end-point  quadrature  formula,  this  singularity  is  peculiar  to  scheme 
I of  the  method  of  integral  relations  and  cannot  be  removed  by 
utilizing  multi-strip  formulations,  although  the  particular  Mach 
number  at  which  the  singularity  occurs  might  be  different  from  that 
of  the  one-strip  formulation.  If,  on  the  other  hand,  the  governing 
ordinary  differential  equations  are  integrated  again  in  the  r- 
direction,  the  singularity  disappears  since  we  now  have  algebraic 
equations.  This  is  the  scheme  III  of  the  method  of  integral  relations, 
which  will  be  discussed  after  we  complete  our  consideration  of  the 
wall- jet  region  in  the  scheme  I formulation. 


B.  V'ALL-JET  REGION.  Integrating  Equations  (2)  and  (7)  from 
the  place  to  the  upper  boundary  of  the  wall  jet,  we  obtain 


^[r^ep.u.v.]  + 2r^8(p. 
dr‘  j 1 1 3 

and 

a ( 1 r 2 

-9-) 

Because  of  Equation  (14) , these  governing  equations  are  considerably 

simpler  than  the  corresponding  ones  in  the  shock  layer.  Utilizing 

Equations  (15),  (16)  and  (30),  one  can  conclude  that  Equations  (14), 

(31)  and  (32)  are  the  governing  equations  for  the  variedales  e,  6 

and  u . Initial  conditions  are,  at  r = 1 
w 


Pw^  = 


(31) 


1/(Y-1) 


li- 


(32) 


12  I 

! 

I 


*7^ 
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e = e. 


a = u , 
w wl 


6 = 6. 


The  first  two  are  supplied  by  the  shock-layer  solution/  and  the  third 
by  using  Equation  (9)  and  the  shock-layer  solution. 


Note  that 


d_ 

dr 


r , l/(T-in  , (2-r)/(Y-l)r  /^  + l\  2l 

[u„(l  - u2)  J = <1  - ul)  [l  . (^)ujj 


du 


w 


dr 


Equation  (32)  becomes  singular  as 


“w  = (^) 


1/2 


= u. 


wn 


(33) 


Utilizing  the  energy  equation  and  the  definition  of  the  speed  of 
sound 


(33a) 


one  may  show  that  Equation  (33)  implies  that 


u = a 
w w 


(33b) 


Therefore,  the  singular  point  is  the  surface  sonic  point,  r = n. 
Since  the  wall  velocity  at  r = n is  continuous  for  a smooth  plate, 
we  may  impose  the  regularity  condition  that,  at  r = n 

V(Y-l)  n 


1 + csc6^(l  - cot  6^) 


u_(l  - u‘  ) 


wn 


wn 


qjd  “ qj) 


T7TFTT 


(- 


- cot6  ) 

n 


2 ’ Pwn^ 

23cot6^csc^6^ — J 


= 0 


(34) 


p4q 


j'^j 


so  that  gj—  is  finite  there.  The  subscript  n denotes  quantities 

evaluated  at  r =»  n.  Equation  (34),  derived  after  some  tedious  but 
straightforward  algebra  from  Equations  (31)  and  (32) , provides  the 
missing  initial  condition  of  the  shock-layer  equations.  This 
completes  the  formulation  of  the  scheme  I of  the  method  of  integral 
relations . 
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METHOD  OF  INTEGRAL  RELATIONS  - SCHEME  III 

Since  Gununer  and  Hunt  TRef"!  (8) ) could  not  'rind  solutions,  that 
will  satisfy  the  wall-jet  equations  by  the  scheme  I of  the  method  of 
integral  relations,  and  since  they  and  South  (Ref.  (11))  have  pointed 
out  the  singular  behavior  of  Equation  (28)  for  low  supersonic  Mach 
numbers,  the  scheme  III  of  the  method  of  integral  relations  is  used 
in  the  present  study.  Two  different  formulations  have  been 
considered  and  they  will  be  discussed  in  the  following. 

A.  ONE-BY-TWO  SOLUTION . Consider  first  the  simplest  case  that 
the  flow  field  between  r ■=  0 and  r = n is  divided  into  two  zones: 

0 ^ r i 1 and  1 ^ r 5 n.  Consider,  in  the  shock  layer,  the 
simplest  approximation 


= CO  to  " r cotOj^ 


which  can  be  integrated  to  yield 


2 

r cot0. 


e = e, 


where  = e (r  = 0),  Equation  (36)  gives  the  relation  between  the 


shock  distances  and  as 


coto^  = 2(Eq  - Sj^)  (37) 

Integrating  Equation  (18)  from  r = 0 to  1 and  utilizing  Equation  (19) , 
we  obtain 

1 

/ r^{3(Pg  - p^)  + PgVg^'^s  UgCota))dr  = 0 (38) 

0 

The  terms  inside  the  curly  brackets  are  even  functions  of  r.  Hence, 
we  may  use  the  simplest  approximating  fxinction 

f(r)  = + (f^  - f(,)r2 


and  Equation  (38)  thus  becomes 


'=sl“sl''slh 


+ 8'PsO  - PwO>] 


+ T3-+-3T  + 6(Pgi  - = 0 (39) 

Obviously,  Equation  (39),  being  an  algebraic  equation,  is  nonsingular. 
Similar  application  of  the  simplest  approximating  function  to 
Equation  (29)  yields 


I ^ 


Vf- 
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2 1/(Y-1) 


^ 


1/(Y-1) 


2 1/(Y-1)  2 V(Y~1) 

so)  2(1  -q2^) 


(j  + 1)  (j  + 3) 


rr+  3) 


= 0 (40) 


We  could  use,  instead  of  Equation  (40),  an  equation  of  global  mass 
conservation 


P V 

00 


Obviously,  is  related  to  u^^^  by  Equations  (4)  and  (30)  as 

r,  „2 


2 1"^ 
1 -<i 

^0 


Note  that  Equation  (41)  is  independent  of  the  approximating  functions 
used  in  the  radial  direction.  It  depends  only  on  the  assumption  of 
a linear  variation  of  pu  with  y,  which  is  always  the  case  for  a one- 
strip  formulation. 


In  the  wall  jet,  1 £ r n»  consider 


- = cot6  ^ ^ — =-  [(1  - r2)cot6„  + n(r2  - n2)cot6,]  (42) 

n(l  - n^)  ^ 

which  yields,  after  a straightforward  integration  process 

2 

e = e,  + [2(n^cot6,  - cot6  ) 

^ 4n(l  - n^)  1 n 

+ (r^  + 1) (cot6^  - ncotS^)]  (43) 

which  gives  the  relation  between  and  6„as 

n n 

2 

— -■y--— )-  (ncot6j^  + cot6^)  (44) 


Equation  (31)  can  be  integrated  from  r * 1 to  n to  yield 
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- (n^£:^sin6^cos6^  - e^sin6^cos6 +23^^ (p ^ - P^)dr  = 


0 (45) 


Consider  the  simplest  approximation  that 

= [(n^  - r2)p^^  + (r^  - l)P^^J/(n^  - D 
Equation  (45)  thus  becomes 

- p^q^  (n^ej^sin6^cos6^  - £j^sin5j^cos6 ^^)  + 23|pjkj^ 
IP„l!n\  - kj)  + P„„(k2  - kj)]  . 

(n^  - 1)  ! 


where 


- ll/(j  + 1) 


^2  “ - 13/(j  + 3) 


Similarly,  Equation  (32)  yields 


2 1/(Y-1) 

q^ci  - qp 


/y  - l\^/2/  2 \V(Y-1)1 

* (M)  (ffr)  J 


= ejq.d  - q.) 


2 1/(Y-1) 


2 1/(Y“1)' 

sin6i  4.  u^,^(l  - u^^) 


The  basic  governing  nonlinear  algebraic  equations  for  the  one-by- 
tv’o  formulation  are  Equations  (39),  (40)  or  (41),  (46),  (47)  and 
(34)  for  the  five  basic  unknowns;  Cq,  n,  and  6^.  We  note 

that  it  is  the  consideration  of  the  surface  sonic  point  which 
provides  two  conditions  (Eqs,  (33)  and  (34)  at  r = n)  with  one 
unknown  (the  location  of  n)  that  enables  us  to  close  the  system. 

We  shall  designate  solutions  obtained  from  using  Equation  (40) , the 
modified  continuity  equation,  by  the  symbol  MCE,  and  those  from 
Equation  (41) , the  global  mass  conservation  equation,  by  the  symbol 
CMC. 

B.  ONE-BY-TH^E  SOLUTION.  In  tliis  formulation  the  wall- jet 
region  is  not  modified.  The shock  layer  is  divided  into  two  regions ; 
1 1 

0 - r i j and  ^ r ^ 1.  Denote  the  quantities  evaluated  at 
r *»  j by  the  subscript  2 and  consider  a continuous  approximating 
function 
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- ^ = cota  ~ r[8(l  - r^)coto^  + (4r^  - l)coto^j/3  (41 

Direct  integration  yields  the  equation  of  shock  detachment  distance 

2 2 
£ = - r I (8cot02  “ cotCj^)  + 2(cot0j  ~ 2cota2>r  ]/6 

After  some  algebra,  one  may  obtain  the  following  relations  between 
the  shock  angles  and  the  detachment  distances: 


cota,^  = (9e0  “ “ 8e2)/3 

cotOj  = (32e,  - 14Ej^  - 18£q)/3 
Equations  (28)  and  (29)  are  of  the  form 

g + r^g  = 0 

where  g is  an  even  function  of  r.  Therefore,  one  may  obtain  by 
straightforward  integrations  that  ^^-2 


^2  ~ ^0 


f - f + 
^1  ^0 


gdr  = 0 


gdr  =■■  0 


The  even  function  g may  be  approximated  by  the  Lagrangian  inter- 
polation formula 

g ' g^d  - r^)  (1  - 4r^)  + gj^(4r^  - l)r^/2  + ISgjd  - r^)r^/3  (54) 

so  that  the  integrals  in  Equations  (52)  and  (53)  become 

. 

J r^gdr  ==  (H^gp  + + H2g2)  (52a) 


gdr  = Iq^o  ^l^l  + ^2^2 


where 


B - i-.  5 . 1 

"o  "*  TT+  ll  4 (j  + ST  413"^) 


«1  = 


- 
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H - (3i  ->■  17) 

”2  " TH  + 3}  (d  + 5) 


^0  - “ TTT~JT  (J  T sT 


T = (3i  -f  7) 

■1  3(j  + 3)0  + 5) 


and 


I.  = 


32 


2 0 + 3)0 '+  sr 


(52d) 

(53b) 

(53c) 

(53d) 


We  tlierefore  have  four  nonlinear  algebraic  equations  obtainable  from 
Equations  (28)  and  (29) . In  addition,  there  are  Equations  (46) , (47) 
and  (34)  of  the  wall-jet  region.  We  now  have  two  additional  basic 
variables,  namely,  e,  and  u^j*  system  is  again  closed.  This 

formulation  is  termed  the  one-by-three  MCE  method.  0ns  may  also 
consider  a one-by-three  GMC  method  by  using  Equation  (41)  to  replace 
the  equation  obtained  by  integrating  Equation  (29)  from  r = 0 to  1. 

It  is  obvious  that  other  approximating  functions  can  also  be 
used.  For  example,  if,  instead  of  the  continuous  representation  as 
given  by  Equation  (54),  the  even  function  g is  assumed  to  be  only 
piecewise  smooth  such  as 

g s gQ  + 4r^(g2  - g^)  for  0 < r < 


and 


g ~ THg,  - g,  + 4r^(g,  - g-)3  for  i < r ^ 1 


Equations  (J2)  and  (53)  still  hold  but  the  constant  coefficients,  H's 
and  I's,  will  be  modified  accordingly.  This  constitutes  the  one-by- 
three  MCE-PWS  method  and  the  corresponding  one-by~ three  GMC-PWS 
method.  Of  course  Equations  (48)  to  (50)  will  also  be  replaced  by 
the  following  piecewise  smooth  equations; 


For 


0 < r < i 


- r cot62 


emd  for  j S r ii  1 


e « £2  - (2(8  cotfij  - cot6j^)  + (4r^  + 1)  (cotf^  - 2 00^2)  3 


where 


cot6. 


4(£o  - 
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cot6^  = j(5e2  "■  - ^^0^ 

Different  approximating  functions  can  also  be  used  in  the  one-by-two 
method.  One  possible  utilization  is  illustrated  in  the  following 
consideration  of  the  stagnation-point  quantities. 

C.  STAGNATION-POINT  VELOCITY  CLIENT.  Of  particular  interest 
to  us  is  the  stagnation-point  velocity  gradient  which  is  directly 
related  to  the  heat-transfer  rate.  Since  u^  is  determined  only  at 

discrete  locations  in  the  scheme  III  of  the  method  of  integral 
relations,  differentiation  of  an  interpolation  formula  is  not 
accurate.  This  difficulty  can  be  circumvented  by  the  following 
method. 

Dividing  Equation  (28)  by  r^  and  taking  the  limit  as  r 0, 
we  obtain 

^ {®^PsO  “ ^wO^  ^sO'^so}  ° 
Similarly,  Equation  (29)  yields 


2 1/(Y-1) 


(&l  • (&), 


2 1/(Y-1) 


Eliminating  (du  /dr)  from  the  above  two  equations,  we  obtain 
® 0 


(S').  ■ 


26(1  - v^g) 


2 1/(Y-1) 


^PsO  ■'  PwO^ 


(1  + 3)P 


sO'^sO^O 


At  r = 0,  a = IT/2.  From  Equations  (23)  to  (26),  (30)  and  (55),  one 
may  conclude  that,  for  fixed  values  of  and  y,  the  stagnation- 
point  velocity  gradient  is  inversely  proportional  to  the  shock 
detachment  distance  at  the  stagnation  point.  Figure  2 shows  the 


value  of  (1  + j) 


as  a function  of  for  y = 1.4, 


Since 


W /o  (Y  + Wo 


we  may  also  obtain  the  relation  that 
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°o 


(Y  ■>■  l)Mf[8(P30  - P„o)  -f  P^qV^;,! 

(1  + - l)P3o''so''»^0 


(56) 


Equation  (56)  may  be  used  to  generate  slightly  more  complicated 
equations  for  the  shock-layer  thickness  and  the  shock  angle.  For 
example,  for  the  one-by-two  method,  we  may  replace  Equation  (35)  by 
the  following  more  complicated  function 


- ^ = cota  « r(  (r^  - DcTq  + r^cota^^] 
Equations  (36)  and  (37)  are  thus  replaced  by,  respectively. 


and 


e = Eq  + r^[0Q(2  - r^)  - r^coto^] /4 


cotCj^  = Oq  + 4(Eq  - e^) 


(57) 

(58) 
(53) 


The  forms  of  other  equations  are  unmodified.  This  formulation  is 
termed  the  one-by-two  GMC  (or  MCE)-SP  method.  In  essence,  the  utili- 
zation of  Equation  (56)  has  increased  the  order  of  the  function  by  2. 
For  example.  Equations  (35)  and  (36)  are,  respectively,  linear  and 
quadratic  in  r,  but  Equations  (57)  and  (58)  are  cubic  and  quartic  in 
r,  respectively.  All  one-by-three  methods  can  be  similarly  modified 
by  incorporating  Equation  (56)  in  their  representation  of  the  shock 
angle  and  the  shock  detachment  distance,  and  will  be  termed 
accordingly. 

RESULTS  AND  DISCUSSION 


The  governing  coupled  nonlinear  algebraic  equations  are  solved 
iteratively  by  the  Newton-Raphson  method.  All  of  the  one-strip 
solutions  obtained  so  far  are  tabulated  in  Tables  1 to  3.  Most  oi 
the  results  do  not  go  above  M„  = 4.  This  is  because,  for  shock- 
interference  problems,  we  are  mostly  interested  in  lower  supersonic 
Mach  numbers.  There  is,  however,  an  upper  limit  on  the  free-stream 
Mach  number  above  which  no  physically  acceptable  solutions  can  bs 
obtained  by  the  present  one-strip  formulation  of  the  method  of 
integral  relations.  This  happens  when  the  location  of  the  surface 
sonic  point,  n#  is  along  the  line  of  the  jet  edge  (r  - 1) . The  trend, 
that  n decreases  toward  unity  as  M«  increases  as  predicted  by  the 
theory,  was  also  observed  experimentally  by  Hunt  and  co-workers 
(Refs.  (8)  and  (14)).  However,  the  actual  occurrence  of  n = 1 io 
believed  to  be  due  to  the  approximation  introduced  by  the  solution 
method.  Fortunately,  this  generally  occurs  above  M<»  = 4 and  hence 
is  not  of  serious  concern  to  us  for  the  present  problem. 


(14V  Carling,  J.  C.  and  Hunt,  B,  L.,  "The  Near  Wall  Jet  of  a 

Normally  Impinging,  Uniform,  Axisymmetric,  Supersonic  Jet," 
Journal  of  Fluid  Mechanics,  Vol.  66,  1974,  pp.  159-176 
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There  is  also  a lower  lisnit  on  M*  below  which  no  physically 
acceptable  solutions  can  be  obtained.  For  the  planar  case,  this 
happens  when  the  calculated  value  of  reaches  unity.  The  fact 

that  it  occurs  at  M„  > 1 Is  again  due  to  the  approximate  nature  of 
the  solution  method.  For  the  axisymmetric  case,  this  happens  at  a 
much  higher  value  of  M«,/  and  the  reason  for  its  occurrence  is  not 
understood  at  the  present  time.  Fortunately,  a quite  wide  range  of 
Moo  does  2xist  between  which  meaningful  solutions  have  been  obtained. 
Because  of  this  much  higher  value  of  the  lower  limit  on  Moo  for  the 
axisymmetric  case,  the  majority  of  the  results  obtained  is  for  the 
planar  case  and  these  results  will  be  discussed  first.  The  results 
for  axisymmetric  flows  will  be  briefly  considered  later.  All  results 
shown  are  for  y - 1.4. 


PLANAR  JET  IMPINGEl^NT 

The  results  of  the  stagnation-point  velocity  gradient  as  obtained 
by  the  various  methods  are  shown  in  Figures  3 and  4 as  a function  of 
Moo-  All  solutions  show  the  same  trend,  namely,  the  initial  rapid 


increase  of 


at  low  Mach  numbers. 


and  the  slow  rise  toward  the 


asymptote  at  high  Mach  numbers.  The  difference  between  onc-by-two 
and  one-by- three  formulations  is  seen  to  be  moderate  at  high  Mach 
numbers,  and  it  drops  very  rapidly  as  M^  is  decreasec'.  The  same  can 
be  said  in  regard  to  the  different  choice  of  the  governing  equations 
between  GMC  and  MCE  methods.  The  application  of  more  complicated 
profiles  (SP  method)  greatly  reduces  the  differences  between  one-by~ 
two  and  one-by-three  formulations,  but  one-by-three  results  display 
only  small  effects  by  the  application  of  these  more  complicated 
profiles,  in  fact,  results  indicate  that  the  one-by- three  formu- 
lation is  quite  insensitive  to  different  approximating  functions 
employed  in  general.  This  is  not  always  the  case  when  other 
quantities  away  from  the  stagnation  point  are  considered,  as  we  shall 
see  later. 


The  detachment  distance  of  the  shock  and  the  upper  boundary  of 
the  wall  jet  as  predicted  by  the  corresponding  one-by-two  and  one-by- 
three  formulations  is  shown  in  Figures  5 to  7 according  to  different 
applications  of  the  method  of  integral  relations.  All  results  show 
the  following  trend;  (1)  both  the  shock  layer  and  the  wall- jet  layer 
become  thicker  as  Moj  decreases;  (2)  aa  Mm  decreases,  the  location  of 
the  surface  sonic  point  moves  away  from  the  line  of  the  jet  edge 
(r  = 1) ; and  (3)  for  a fixed  Mm,  the  moderate  difference  between 
one-by-two  and  one-by-three  formulations  at  the  symmetry  line 
{r  = 0)  is  reduced  even  further  at  the  line  of  the  jet  edge  (r  = 1) . 

The  surface  Mach  number  evaluated  at  r = 3. , and  the  Mach 

number  behind  the  shock  at  r = 1,  are  depicted  in  Figure  8 as 

functions  of  M^.  Clearly,  neither  nor  M^j^  is  generally  equal  to 

unity.  Hence  the  boundary  conditions  employed  in  References  (8)  and 
(9)  are  incorrect.  The  corresponding  values  of  the  shock  angle  at 
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the  line  of  the  jet  edge,  as  obtained  from  various  methods  are 
shown  in  Figures  9 and  10.  Similar  to  they  are  seen  to  be  more 

method-dependent  than  quantities  such  as 

The  surface  pressure  distribution,  as  shown  in  Figure  11,  indi- 
cates the  general  insensitivity  of  the  results  to  various  schemes 
employed.  The  only  noticeable  difference  is  the  somewhat  fuller 
profile  predicted  by  the  one-by- three  formulation. 

It  therefore  appears  from  self-consistency  that  reasonable 
engineering  solutions  for  the  stagnation-point  velocity  gradient 
(hence  Eq)  and  (hence  u^^  and  p^j^)  have  been  obtained.  Since 

heat-transfer  rate  is  proportional  to  the  square  root  of  the  velocity 
gradient  at  the  stagnation  point  (Refs.  (15)  and  (16)),  peak-heating 
prediction  is  thus  even  less  method-dependent.  This,  however,  is  in 
direct  contras v.  to  the  axisymmetric  case  which,  to  be  discussed 
next,  is  seen  to  be  far  from  converging. 


AXISYMMETRIC  JET  IMPINGEMENT 

Among  all  the  methods  employed , only  one-by-two  CMC  and  MCE 
schemes  have  produced  solutions  which  appear  not  to  violate  some  of 


the  obvious  physical  constraints  such  as  p - > p 


and. 


wO  *^w2  '^wl  ^wn 
as  Mm  decreases,  both  (du^dr)^  and  u^^^  will  also  decrease.  The 

results  are  tabulated  in  Tables  3a  and  3b.  The  lowest  Mo>  shown  in 
each  table  is  the  lower  limit  of  the  Mach  number  below  which  no 
solution  is  obtainable.  As  we  can  see,  the  corresponding  M^^^  is  far 

from  being  unity.  The  reason  for  the  existence  of  this  relatively 
high  value  of  the  lower  limit  of  M^  is  not  understood  at  the  present 
time. 


The  axisymmetric  results  are  qualitatively  similar  to  the  planar 
solutions.  There  are  noticeable  differences  also.  For  example,  for 
the  axisymmetric  case,  the  shock-layer  thickness  drops  off  at  a much 
faster  rate  as  one  moves  away  from  the  stagnation  point.  This  results 
in  a smaller  shock  angle,  and  a thinner  wall- jet  layer.  In  fact, 

the  rate  that  drops  with  respect  to  decreasing  M^  is  so  large  that 

Mgj^  turns  out  to  be  increasing  slightly  as  M^^  is  decreased.  This 

trend  is  clearly  opposite  to  that  of  the  planar  case  which  shows  the 

(15)  Cohen,  C,  B.  and  Reshotko,  E.,  "Similar  Solutions  for  the 
Compressible  Laminar  Boundary  Layer  with  Heat  Transfer  and 
Pressure  Gradient,"  NACA  Rpt  1293,  Lewis  Research  Center, 
Cleveland,  Ohio,  1956 

(16)  Fay,  J.  A.  and  Riddell,  F.  R. , "Theory  of  Stagnation  Point  Heat 
Transfer  in  Dissociated  Air,"  Journal  of  the  Aeronautical 

Sc' «^nces,  Vol.  25,  1958,  pp.  73-85,  l2l 

22 


NSWC/WOL/TR  75-195 


monotonic  decreasing  behavior  as  was  depicted  in  Figure  8.  Since  the 
dxisymmetric  solution  appears  to  be  very  method-dependent  (as  can  be 
seen  easily  by  the  fact  that  even  one-by-two  CMC-  and  MCE-SP  methods 
yield  no  physically  acceptable  solutions) , results  obtained  by  other 
methods  are  needed  before  these  different  trends  can  be  ascertained 
or  refuted. 

CONCLUSIONS 

The  major  conclusion  that  we  may  draw  from  the  present  study  is 
that  solutions  that  satisfy  all  well-posed  boundary  conditions  can 
be  obtained  by  the  one-strip  formulation  of  the  method  of  integral 
relations.  The  application  of  the  scheme  III  of  the  method  has 
enabled  us  to  avoid  both  the  unwanted  singularity  at  the  low  super- 
sonic Mach  number  and  the  numerical  difficulty  of  satisfying  the 
regularity  condition  at  the  surface  sonic  point  peculiar  to  the 
scheme  I of  the  method.  Rational  engineering  solutions  for  the 
stagnation-point  velocity  gradient  and,  hence,  the  peak  heat- transfer 
rate  have  been  obtained  for  a planar  supersonic  balanced  jet 
impinging  normally  on  a flat  surface.  However,  more  theoretical 
and/or  experimental  studies  are  needed  before  present  results  can  be 
quantitatively  assessed.  Toward  this  goal,  a two-strip  formulation 
of  the  method  of  integral  relations  has  been  completed.  Unfor- 
tunately, because  of  the  time  limitations,  no  quantitative  results 
have  yet  been  obtained.  For  the  sake  of  completeness,  this  formu- 
lation is  included  in  the  Appendix. 

Since,  for  impingement  angles  betv;een  normal  (90  degrees)  and 
about  50  degrees,  the  effect  of  the  angle  of  impingement  on  the 
peak  pressure  was  found  experimentally  by  Henderson  (Ref.  (17))  to  be 
small,  the  present  planar  jet-impingement  model  might  be  coupled  with 
the  shock-interference  model  of  Edney  (Ref.  (1))  as  programmed  by 
Morris  and  Keyes  (Ref.  (18))  to  predict  type  IV  shock-interaction 
effects.  In  view  of  the  extremely  short  computer  time  required  by 
the  present  method  (typically  less  than  five  seconds  on  a CDC  6500 
computer  for  one  converged  solution  at  one  Mach  nximber) , this 
approach  is  indeed  very  attractive. 


Tl7)  Henderson,  L.  P.,  "Experiments  on  the  Impingement  of  a Super- 
sonic Jet  on  a Flat  Plate,"  ZAMP,  Vol.  17,  1966,  pp.  553-569 

(18)  Morris,  D.  J.  and  Keyes,  J.  W. , "Computer  Programs  for  Predict- 
ing Supersonic  and  Hypersonic  Interference  Plow  Fields  and 
Heating,"  NASA  TM  X-2723,  Langley  Research  Center,  Hampton,  Va. 
1973 
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Table  1 PLANAR  JET  IMPINGEMENT:  ONE-BY-THO  SOLUTIONS 

a.  GMC  Method 
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APPENDIX  A 


TWO-STRIP  FORMULATION  OF  THE  JET- IMPINGEMENT  PROBLEM 


METHOD  OF  INTEG^L  RELATpNS  - SCHEME  I 

The  flow  field  is  divided  into  two  strips  in  the  axial  (y-) 
direction  by  the  middle  line  y = e/2.  The  governing  equations  are 
different  depending  on  whether  they  are  in  the  shock  layer 
(0  ^ r i 1)  or  in  the  wall-jet  layer  (1  i r ^ n) . 

A.  SHOCK-LAYEP.  REGION.  Integrating  the  axial  moment  equation 
(2)  from  0 to  e/2,  we  obtain 


(r)/2 


r^puvdy  - I If 


Ph'^H 


- p v2\  = 
^w  w f 


(Al) 


where  the  subscript  H denotes  quantities  evaluated  at  y = e/2.  If  a 
quadratic  profile  in  y is  assumed  for  the  integrands  in  Equations 
(Al)  and  (27),  after  some  algebra,  we  may  obtain,  as  an  approximation 
to  Equation  (2),  the  following  two  ordinary  differential  equations: 


|j-(r3cp^u^v^)  + 4r3(P3V^(v^  + u^coto) 

■ Vh‘^''h  + Pw>)  = “ 


Equations 

Similarly, 


®F<p^"PhPh''k>  * f-  •*  ^ 

2PjjVjj(2Vy  + UjjCOtu)  + 8(Pg  + 4pjj  - 5p^)  } = 0 

(5)  and  (17)  have  been  used  in  the  above  equations. 
Equations  (i)  and  (3)  yield 

a?  [p^‘^<pw“w  - Ps“s>  ’Ph'2''h  * V°“> 

- p (v  + u„coto)}  = 0 
s s s 


^ j^r3  c(2p„Ug  + PgUg)  + r3{5Pg(v-3  + u^coto) 


- 2pjj(2Vjj  + Uj^cotc)}  = 0 


(A2) 


(A3) 


(A4) 


(A5) 


A-1 
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+ ®P«  - ‘’s“s  - 8Ps>]  " + “b““> 

- PgUgCVg  + UgCoto)  + 6(Pjj  - Pg)coto}  = jBe(p^  - p^)  (A6) 

and  ^ 

gjr  j^r^e(PsU^  + 6p^  + 2p„u^  + 28PH>j  + r3{5p^u^(v^  + u^ootp) 

- 2PjjUjj(2Vjj  + UjjCOto)  + 6(5pg  - 2pjj)cota}  = j6e(2pjj  + p^)  <A7) 

The  energy  Equation  (4)  gives 

Pw  = P«<1  - “w> 

and 

Pr  ■ 

Thus,  Equation  (A2)  defines  the  rate  of  changfe  of  o,  Equations  (A4) , 
(A6)  and  (A8)  thore  of  u , p and  p , and  Equations  (A3) , (A5) , (A7) 

Mi  Tt  Mi 

and  (A9)  those  of  Ujj,  and  Pj^.  Because  of  Equation  (18),  one 

may  replace  Equations  (A6)  and  (A8)  by  the  simpler  algebraic 
relations.  Equation  (30)  and 

r, 


1 - u" 


(AlO) 


L J 

Therefore,  there  are  six  ordinary  differential  equations  (Eqs.  (5) , 
(A2)  to  (A5)  and  (A7))  for  a,  e,  u^,  Ujj,  Vjj  and  Pj^.  Initial 

conditions  are,  from  Equations  (19)  to  (21) , at  r = 0 


"ho  ' “ 


"wo  * " 


2 -|1/(Y-1) 


Pbo=  ~ 

and  = tt/2 

The  two  missing  initial  conditions  for  and  v^q  are  supplied  by  the 

two  regularity  conditions  at  the  surface  sonic  point  and  the  singu- 
larity on  the  middle  line.  It  is  known  for  the  jet-impingement 
probl^  that  the  surface  sonic  point,  r « n»  is  outside  the  shock 
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(8)),  but  the  location  of  the  other  "sonic"  point  on  the 
middle  line  relative  to  the  line  of  the  jet  edge  (r  = 1)  is  unknown 
a priori.  These  singularities  and  the  associated  regularity 
conditions,  well-known  for  blunt-body  problems,  are  essential  for 

closing  the  system  of  equations.  They  will  be  discussed  in  detail 
later . 


the  structure  of  Equation  (A2)  is  similar  to  Equation 
(28) , Equation  (A2)  will  also  become  singular  as 


a(PgU^v^) 


do 


= 0 


and  do/dr  will  become  unbounded.  A formulation  based  on  scheme  III 
will  thus  be  required.  Before  this,  however,  we  shall  complete  the 
present  discussion  of  the  scheme  I method  by  considerina  the  wall- 
let  region. 


. ^GION . Integrating  Equations  (1)  to  (3)  from  the 

piate  to  the  middle  line  and  from  the  plate  to  the  upper  boundary 
of  the  wall  jet,  anu  after  some  straightforward  algebra,  we  obtain 


^j,(r^epjUjVj)  + 4r^{6(pj  - + p^)  - PjjVjj(2Vg  + UjjCot6)}=0  (All) 


dr^^^^PnVH^  + |“^2pjjVjj(2Vg  + UgCOt6)  + B(p^  + 4p^  - 5p^) } = 0 (A12) 


gj:[r3e(p^u^  - p^Uj)]  + 4r^py(2vjj  + UjjCotS)  = 0 


(A13) 


[r^e  (2PjjUjj  + Pj^jJ^  " = 0 


(A14) 


d 

dr 


j^r^e  (PjU? 


2 2 1 ’ 
j ®Pj  20Pjj)  + r^{S(5pj  - 2pjj)cot6 


- 2PjjUjj(2Vjj  + UjjCOt6)}  = j0e(2Pjj  + p.) 


(A15) 


As  in  the  shock  layer,  the  other  ordinary  differential 
equation  that  comes  from  the  radial  momentum  equation  (3)  is  re- 
placed by  the  algebraic  equations  (30)  and  (Aid) . In  addition, 
there  is  the  geometric  relation#  Equation  (6) # the  boundary  conditions 
at  the  wall- jet  boundary,  Equations  (8)  to  (16) , and  the  energy 
equation  (A9) . Therefore,  there  are  six  ordinary  differential 

equations  (Eqs.  (6),  and  (All)  to  (A15) ) for  e,  6,  u , u„,  v„  and 

W H H 


A- 3 
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p^.  Matching  conditions  at  r = 1 supply  the  initial  conditions. 

One  may  combine  Equations  (A13)  and  (A14)  to  give 

r3c(p^u„  + 4PgUg  + PjU^)  = + 4p„^u„j  + p.u^^)  (A16) 

which,  being  an  algebraic  relation,  can  be  used  to  replace,  e.g.. 
Equation  (A14) . 

C.  REGULARITY  CONDITIONS.  Utilizing  Equation  (AlO)  and  after 
some  straightforward  algebra,  we  may  rewrite  Equation  (A13)  in  the 
form 

du  N, 

aF"  ■ 

where 

' [i  - <H-r>“w] 

To  have  a finite  value  of  du^/dr  at  the  singularity  given  by 

Equation  (33),  we  require  that  -►  0 as  0 at  r = n.  This 

provides  us  with  the  regularity  condition  which,  using  Equation  (All) 
at  r = n to  get  rid  of  d6^/dr  and  after  some  straightforward  algebra, 

becomes 

q.(?^  - (p^^u^^cos2«^  + p^g.sin^S^) 

- 46cos6^(Pj  - 2pjj^  + p^^)  = 0 ' (A17) 

The  location  of  the  singularity  on  the  middle  line  is  again 
unknown  a priori.  Two  different  formulations  are  needed  depending 
on  whether  it  is  larger  than  1 or  otherwise.  Let's  consider  the 
first  case  (henceforth  referred  to  as  Case  W) , and  denote  the 
singularity  to  be  at  r = C > 1* 

From  Equations  (A9) , (A12) , (A14)  and  (A15)  we  may  obtain 

3r  D2 

where 

- (Y  + l)u^  + (Y  - 1) (v^  - 1) 


A- 4 


i 


•V 
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Therefore,  as  D2  0 at  r = we  need  to  impose  the  regularity 

condition  that  N2  -*•  0 at  r = C*  Using  Equation  (All)  at  r = C to  get 

rid  of  d6^/dr  (where  the  subscript  C denotes  quantities  evaluated 

at  C)  and  after  some  tedious  but  straightforward  algebra,  we  may 
obtain  the  regularity  condition 


2 2 

p^q^sin  6^((y  - l)sin6^  - yqj 


. V ■ 


cot6 , 


+ + u„jCOt6j)[c^(l  + 


'hc  hj 


i“hc 


PHjlcotSj  - -r  PhC  614Cg(p.  - 2pj,j 


♦ P«C>  * VHE<Pj  * ^Ph5  ■ “ “ 


(A18) 


where 


= (Y  - l)qjCOs26^ 

Cg  = (1  - Y + 2YqjUg^sin6^)cos6^ 


and  , at  r = ? 


[(^  + 0 (^  ■ ''»f| 


(A19) 


It  is  easy  to  show  that  Equation  (A19)  is  equivalent  to  Ug^  = a^^. 

Let  us  now  consider  the  case  when  the  singularity  on  the  middle 
line  occurs  in  the  shock  layer.  Henceforth,  we  shall  refer  to  this 
case  as  Case  S and  denote  the  singularity  to  be  at  r = ? < 1. 

From  Equations  (A3) , (A5) , (A7)  and  (A9)  we  may  obtain 


where  again 


(Y  + l)u^  + (Y  - 1)  (v^ 


This  is  to  be  expected  since  the  structure  of  the  governing 
equations  in  both  layers  is  similar.  Therefore,  as  0 at 

r = c,  Nj  -*■  0.  Using  Equation  (A2)  at  r = C to  get  rid  of  da^/dr 

(where  the  subscript  c denotes  quantities  evaluated  at  r » C)  and 
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after  some  tedious  but  straightforward  algebra,  we  obtain,  at 
r = C 

1/2 


(H4)(i  ■ ’hc)] 

j[“H<®Ps  * - «»s“s]  - 

+ Pg(Vg  + UgCOto)  [5UgUjj  + 26(VgVjj  - 5)3  + 2Bpjj(2Vjj  + u^cota) 

(1  + 3Vg  - Up)  + 8(5Pg  - 2pjj)UjjCoto  + 2B^Vjj(Pg  + 4pjj  - 5p^) 


(A20) 


coto) 


- j%€(2p„  + p^)/r 


j-[eu 


dp 
H 3o 


•H 

d(p  up  d(p 

® t u..  - 23 


H 3F 


V] 


- <=<*“>  + ^■''s''s'''8  * »s<=°«)  - i>h''h'2''h 

+ UyCoto)  + B(Pg  - 2P|,  + p„)l}  * 0 

where  Equation  (A21)  is  evaluated  at  r » C. 

D.  STAGNATION-POINT  VELCX^ITY  GRADIENT.  Dividing  Equation 
(A2)  by  r^  and  taking  the  limit  as  r 0,  we  obtain 

du  / 

L\  ^ Ain  r — -i. 

'sO  ‘^HO  10 


(A21) 


+ ®<Pso  - =*Pho  * Pwo*}  - ° 

Similarly,  Equation  (A4)  yields 

du . du 


'HO  " ‘^sO^sO^  * ° 


(1  + :j)^ojpwo^anr^Q " Pso^a^r^Qj  ^^^p^qv^ 

Eliminating  (du^/dr)  q from  the  above  two  equations  we  obtain 

4I2PH„Vgo(VHO  - ’.o’  - «<PsO  - ^PhO  + Pwo” 

\ir)  M— vv-  (A22) 


* J>»wo''.o'o 


for  fixed  values  of  and  Y,  Equation  (A22)  Indicates  that  the 
product  eQ(du^dr.Q  depends  also  on  Vjj^. 


A-6 
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METHOD  OF  INTEGRAL  yiATIONS  - SCHEME  III 

A two-by-four  jfonnulation  will  be  presented  as  an  example  below. 
To  extend  the  formulation  to  two-by-n  with  n > 4 is  straightforward, 
but  the  algebra  involved  will  be  much  more  complicated.  In  addition, 
there  will  be  more  equations  to  solve.  This  certainly  will  aggravate 
the  convergence  problem.  For  simplicity,  we  shall  only  present  the 
details  of  Case  W.  The  other  case  is  very  similar. 

The  flow  field  is  divided  in  the  radial  direction  into  (0,  j, 

1,  5,  n).  In  the  shock  layer.  Equations  (48)  to  «53)  obviously  still 
hold.  In  addition.  Equations  (a2)  to  (A5)  are  of  the  form  of 
Equation  (51) , and  hence  they  can  be  put  into  the  forms  of 
Equations  (52)  and  (53) , with  the  integrals  and  coefficients  given 
by  Equations  (52a)  to  (53d) . Equation  (A7)  is  of  the  form 


dr 


+ r-'g  = jh 


(A23) 


where  g and  h are,  respectively,  odd  and  even  in  r.  In  addition, 
gQ  = 0.  Straightforward  integrations  of  Equation  (A23)  over  r yield 


and 


(A24) 


r-'gdr  = 


(A25) 


Consider  the  continuous  approximating  functions  for  g and  h as 


g = f[892  - + 4(g^^  - 2g2)r^] 

_ (16h2  - h^  - IShg)  . 4{3h^  + h^  - 4h2>  ^ 

h . h(j  + 3 r + 5 r 

Equations  (A24)  and  (A25)  become 

2 


Where 


a, - * a--  “ 0 
10  20 

, . (3j  4) 

*11  3(1  + jjfj  + 4) 

*12  “ 3(1  + J)(l  + '4) 
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®21 

“ 37J 

2"' 

^22 

^ TT 

^10 

b, , 

. 11 

°11 

45 

b,  - 

s 

“12 

45 

b 

» lii 

°20 

60 

^21 

^22 

“ IP- 

j + 


In  the  wall- jet  layer,  consider 


cot6,  + 


S'  y ^ cot6. 

n(5^  - n^)(l  - n^)  5(n^  - r)  (1  - 5^)  ^ 


Direct  integration  yields 


e a e. 


- (y  - 1) 


cot  6. 


2(r^  + r^  + 1)]  + 


2(r^  + r^  + 1)1  + 


(n*'  - 1)  (5‘  - 1) 
cot6^ 

+ 2 

niK  - n^)  (1  - n^) 

cot6- 

+ 2 1 

5(n^  - 5 ) (1  - 


[6n^?^  - 3(n^  + 5^) (r^  + 1)  + 


2(r^  + r^  + Dll 


C65^  - 3(5^  + 1) (r^  + 1)  + 


[6n^  - 3(n^  + 1) (r^  + 1)  + 


(A7& 


i 
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After  some  algebra,  we  obtain 


e = e - 

1 12 


) ( 

((n^  - 1) 


(3n^  - - 2)  - 


(A29) 


cot6. 


n 


- n^)  (1  - n^) 


(r  - 1)  - 


cot  6- 
5 ^ 2~ 

- r) 


(25^  - 3n^  + 1) 


and  e is  obtained  by  interchanging  ^ and  n in  Equation  (A29) . 


Equations  (All)  to  (A14)  are  of  the  form  of  Equation  (51) . The 
even  function  g can  now  be  approximated  again  by  the  Lagrangian 
interpolation  formula 

(n^  - r^)  (C^  - r^)  (g^  - r^)  (1  - xh 

(n^  - 1)  ir  - 1)  ^ (r  - n^)  (1  - n^)  ^ 

. (n^  - r^)  (1  - r^)  _ 

+ — X 5 5“  g-  (A30* 

(n  - C^)(l  - C ) ^ 


which  is  symmetric  in  5 and  n»  i.e.,  the  equation  is  un -.ranged  by 
interchanging  C and  n.  Integrating  Equation  (51)  over  r and  using 
Equation  (A30) , we  obtain 


where 


and 


+ gj^G(n,5,l  ; a)  + g G(l,?,ri  ; a)  + 


g^G(l,T),C  ; a)  = 0 ; a = 5,n 


(A31a,b) 


S. (x,y;r) 

G(x,y,z  ; r)  h — * — * y- 

(x^  - z^)  (y*^  - z^) 


(A3  2) 


Sj (x,y;r) 


2 2 
X y 

Tr+ IT 


(j+i) 


(x^  + yh  r_(i+3) 
• (j  1 1)  [’^ 


1 rr<3«)  . il 

- 

Equation  (A15)  is  of  the  form  of  Equation  (A23) . Using 
Lagrangian  interpolation  formula  for  approximating  the  odd  and  even 
functions  g and  h,  respectively,  we  obtain  by  straightfosrward 
integration 
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f - f.  + g,H(rw5fl  ; a)  + g„H(l,?,n;  a)  + gpH(l,n,C  ; a)  = 
all  n ^ 

j[h,K{ri,^,l  ; a)  + h K(l,^,n  ; a)  + h^Kd^n,^  ; a)l  ; a = 

(A34a,b) 


where 


and 


S.^l(x,y  ; r) 

H(x,y,z  ; r)  = * * 5 x- 

z(x  - z ) (y  - z ) 


K(x,y,z  ; r)  = — j 


Sq(x,y  } r) 


/ **  2^/2  ^ 
(x  - z ) (y  - z ) 


(A3  5) 


(A36) 


There  are  22  basic  unkuowns  in  the  two-by-four  formulation:  e^,  ^2 

El,  6^,  6^,  '^wl'  ^HO'  '^H2"  '^Hl'  ^Hn'  ^H2'  '^Hl' 

^H2'  ^Hn'  ^ basic  equations  are:  Equations 

(A2)  to  (A5)  in  the  forms  of  Equations  (52)  and  (53);  Equation  (A7) 
in  the  forms  of  Equations  (A26a,b);  Equations  (All)  to  (Al4)  in  the 
forms  of  Equations  (A31a,b) ; Equation  (A15)  in  the  form  of 
Equations  (A34a,b);  and  the  regularity  conditions,  Equations  (A17) 
and  (A18).  Total  number  of  equations  is  also  22  and  the  system  is 
closed. 


FLIGHT  MIIASUREKENTS  DIVISION 
EXTERNAL  DISTRIBUTION  LIST  (A-1) 


Copies 

Commander,  Naval  Sea 
Systems:  Connano,  Hqs. 

Department  of  the  Navy 
Nashinqton,  D.  C.  20360 
Chief  Tech.  Analyst 
SEA  05121 
SEA  033 
SEA  031 

SEA  09G32  2 

SEA  035 

Commander,  Naval  Air  Systems 
Command , Hqs . 

Department  of  the  Navy 
Washington,  D.  C.  20360 
AIR  03B 
AIR  03C 
AIR  320 
AIR  320C 

Dr.  H.  J.  Hueller,AIR  310 

AIR  S0174  2 

Office  of  Navy  Research 
800  N.  Quincy  St. 

Arlington,  Va.  22217 
ONR  100 

Horton  Cooper,  430B  2 

Commander 

Naval  Ship  Research  and 
Development  Center 
Bethesda,  Md.  2003S 

Central  Library  Br.  (5641) 

Aerodynamics  Lab.  (5643) 

Coonander,  Naval  Weapons  Center 
China  La)ce,  Calif.  93SSS 
technical  Lib.  (533) 

Code  406 

R.  E.  Meeker  (4063) 

Director,  U.  S.  Naval 
Research  Laboratory 
Washington,  D.  C.  20390 
Library 
Code  6503 

NASA 

Langley  Research  Center 
Hampton,  Va.  23665 

HS/18S  Technical  Library 
Aero  t Space  Kcch.  Div. 

Dennis  Bushr.ell 
Ivan  Beckwith 
R.  Trlmpi 
Julius  Harris 

N^ 

L^is  Research  Center 
21000  Brookpart  Road 
Cleveland,  Ohio  44135 
Library  60-3 

Ch,  Wind  Tunnel  i Flight  Div. 

NASA 

George  C.  ,v.arshaXl  Space 
Flight  Center 
Huntsville,  Ala.  35812 
Mr.  T.  Reed,  R-?ER0-AU 
Mr.  W,  K.  Dahm, 

NASA 

600  Independence  Avo.,  S.  H. 

Washington,  D.  C.  20546 
F.  C.  Schwenk,  Director, 

Research  (Code  RR) 


Copies 

NASA 

P.  O.  Box  33 

College  Park,  Md.  20740 

NASA  Ames  Research 
Moffett  Field,  Ca.  94uj5 
Dr.  M.  Horstman 
P.  Kutlcr 
J.  Rakich 
R.  MacCormack 
L.  H.  Jorgensen 
E.  J.  Hopkins 
B.  H.  Album 
E.  R.  Keener 

Technical  Library 
Director  Defense  Research  and 
Engineering  (DDRtE) 

Room  3E-1063,  The  Pentagon 
Washington,  D.  C.  20301 
Stop  103 

Defense  Documentation  Center 
Cameron  Station 

Alexandria,  Va.  22314  12 

Commander  (5632.2) 

Naval  Missile  Center 
Point  Mugu,  Ca.  93041 
Technical  Library 

Commanding  Officer 
USA  Aberdeen  Research  and 
Development  Certer 
Aberdeen  Proving  Ground, 

Maryland  21005 

8TEAP-TL  (Tech  Lib  Div) 

AKXRD-XSE 


Director,  Strageti,c  Systems 
Project  Office 
Department  of  the  Navy 
Washington,  D.  C.  20390 
SP-2722 

Director  of  Intelligence 
Hdqs.,  USAP  (AFNINDE) 

Washington,  D.  C.  20330 

AFOIN-3B  2 

Los  Angeles  Air  Force  Station 
EAMSO/DYAE 

P.O.  Box  92960,  Worldway  Postal  Center 
Los  Angelos,  CA  90009 
Code  PSSE 
Code  RSSH 

Beadquartcis,  Arnold  Engineering 
Development  Center 

Arnold  Air  Force  Station,  Tenn.  37389 
Library  Documents 
R.  W.  Hensol,  .D 
Capt.  C.  Tirres/DYR 
C.  Welsh 

von  Karman  Gas  Dynamics  Facility 
ARO,  Inc. 

Arnold  Air  Force  Station,  Tann.  37369 
Dr.  J.  Whitfield,  Chief 
L.  H.  Jenke 
W.  8.  Baker,  Jr. 


t 

I 

I 

I 

I 


I 

I 

I 


DISTRIBUTIOK  (CONT'D) 


Conmianding  Officer,  Harry 
Diamond  Laboratories 
Washington,  D.  C.  20438 
Library,  Rm  211,  Bldr . 92 

Connianding  General 
U.  S.  Army  Missile  Command 
Redstone  Arsenal,  Ala.  35809 
amshi-rr 

Ch,  Document  Sec. 

AMSHI-RDK,  Mr.  R.  Deep 
AMSMI-RDK,  Mr.  T.  Street 
D.  J.  Spring 

Department  of  the  Army- 
Office  of  the  Chief  of  Research 
and  Development 
ABMDA,  The  Pentagon 
Washington,  D.  C.  20350 

Commanding  Officer 
Picati.iny  Arsenal 
Dover,  N.  J.  07801 
Nr.  A.  A.  Loeb 
SKOPA-VC-3 

Commander  (ADD 

Naval  Air  Development  Center 

Johnsville,  Pa.  18974 

Air  Force  Weapons  Laboratory 
Xlrtland  Air  Force  Base 
Albuquerque,  N.  M.  87117 
Technical  Library  (SOL) 

Capt.  Tolman/SAS 

0.  8.  Army  Ballistic  Missile 
Defense  Agency 
1300  Wilson  Blvd. 

Arlington,  Va.  22209 
Dr.  S.  Ale.;ander 

The  Johns  Hop)(ins  University 
(C/NOw  7386) 

Applied  Physics  Laboratory 
8621  Georgia  Ave. 

Silver  Spring,  Md.  20910 

Document  Library  2 

Dr.  F.  Hill 
Dr.  L.  Cronvich 

Pii^ector,  Defense  Nuclear  Agency 
Headquarters  DASA 
Washington,  D.  C.  20305 
STSP  (SPAS) 

Commanding  Officer 
Naval  Intelligence  Support  Center 
4301  Suitland  Road 
Washingto.i,  D.  C.  20390 

Department  of  Aeronautics 
DFAN 

USAF  Academy 
Colorado  80840 
Col.  D.  H.  Daley 
Capt.  J.  Williams 

Armament  Development  and  Test 
Center 

Eglin  AFB,  Fla. 

Technical  Lib,  DLOSL 

Headquarters , Edgewood  Arsenal 
Edgewood  Arsc.-.al,  .Md.  21010 
A.  Flataii 


Copies 


2 


Cosnander 

Natic)c  Development  Center 
Natic)c,  Mass.  01760 
AKXNH-UBS 
G.  A.  Barnard 


AFFDL/FX 

Hright-Patterson  Air  Force  Base 
Dayton,  Ohio  45433 
Dr.  D.  J.  Harney 

AFFDL/FXC 

Wright-Patterson  Air  Force  Base 
Dayton,  Ohio  45433 
Mr.  M.  Buck 
P.  Ciragosian 

Naval  Air  Test  Facility 
Lake  Hurst,  N.  j.  08733 
Dr.  H.  Sule 

**my  Aviation  Systems  Command 
p.  O.  Box  209,  Mein  Office 
8t.  Louis,  Mo,  63166 
Dr.  L.  Lijewski 


Copies 


2 


■V 


I 


I 

I 

1 


experimental  aerodynamics  branch 

EXTERNAL  DISTRIBUTION  LIST  (A2) 


I 


i 


i 

i 


CcpieM 

Aerospace  Engineering 
Program 

Univcroity  of  Alabama 
P.  O.  Box  2908 
Univernity,  Alabama  35406 
P*-o£.  W.  X.  Key,  Chm. 


Cople* 

Caae  Wcctern  Reaerve 
Univnrsity 

Division  of  Fluid,  Thermal 
and  Aerospace  Engineering 
Cleveland,  Ohio  44106 
Or.  Eli  Keshotko,  Head 


AME  Department 
University  of  Arizona 
Tucson,  Arizona  85721 
Dr.  L.  D.  Scott 

Polytec^aic  Institute  of 
Brooklyn 

Graduate  Center  Library 
Route  110,  Farminqdale 
Long  Island,  New  York  11735 
Dr.  .1  Polczynski 

Polytechnic  Institute  of 
Brooklyn 
Spicer  Library 
333  Jay  Street 
Ercoklyn,  New  York  11201 
Reference  Department 


California  Institute  of 
Technology 
Pasadena,  C'  91109 
Graduate  -ieronautlcal  Labs. 
Aero.  Librarian 
D.  Coles.  321 
Guggenheim  Lab. 

Or.  A.  Roshko 


jmiversity  of  California 

»«'»anioal  Engineering 
Berkeley,  CA  94720  ’ 

Prof.  R.  Grief 


Motre  Dame  University 
Votre  Dame,  Indiana  46556 
Dr.  V.  Goddard 
Dr.  V.  Nee 
Dr.  T.  Muller 
Ox  R.  Nelson 
Dr.  P.  Ravan 
Di*o£.  R.  .Elkenberry 


CASDYNAMICS 

Oniversfty  of  California 
Richmond  Field  Station 
1301  South  46th  Street 
Ririsuond-  C...l'ornia  94804 
A.  K.  Opponheim 

Department  of  Aerospace 
Engineorlng 
University  of  Southern 
California 
Univercitv  Park 
Los  Angeles 
California  90007 
Dr.  John  I.aufvr 

University  of  California  - 
San  Dieno 

Department  of  Aerospace  and 
Mechanical  Engineering 
Scierces 

LaJolla,  California  S2037 
Dr.  P.  A.  Libby 


The  Catholic  University  of 
Anerics 

Washington,  D.C.  20017 
Dr.  C.  C.  Char.g 
Dr.  Paul  K.  Chang 

Mechanical  Engr.  Dept. 
Dr.  M.  J.  Casarella 
Mechanical  Engr.  Dept. 

University  of  Cincinnati 

Cincinnati,  Ohio  45221 
Department  of  Aerospace 
Engineering 
Dr.  Arrold  Polak 

Department  of  Aerospace 
Engineering  Sciences 

Onlversity  of  Colorado 

Eoulder,  Colorado  80302 

Cornell  University 

Graduate  School  of  Aero. 
Engineering 

Ithaca,  New  York  14850 
Prof  H.  R.  Sears 
Dr.  S.  F.  Shen 
Prof.  P.  K.  Moore,  Head 
Thermal  Engineering 
Dept.,  208  Upson  Hall 

University  of  Delaware 

Mechanical  and  Aeronautical 
Engineering  Dept. 

Dewark,  Delaware  19711 
Dr.  James  £.  Danberg 

Georgia  Institute  of 
Technology 

225  North  Avenue,  N.W. 

Atlanta,  Georgia  30332 
Or.  Arnold  L.  Ducoffe 


Technical  reports 
Collection 
Gordon  McKay  Library 
Harvard  University 
Div.  of  Eng'g.  and 
Applied  Physics 
Pierce  Hall 
Oxford  Street 
Cambridge 

Massachusetts  02138 

Illinois  Institute  of 
Technology 
3300  South  Federal 
Chicago,  Illinois  60616 
Dr.  K.  V.  Morkovin 
Prof.  A.  A.  Fejer 
N.A.E.  Dept. 

University  of  Illinois 
ICl  Transportation  Bldg. 
Urbana,  Illinois  61801 
Aeror.xutical  and 
A>  t.c  inautlcal 
Engineering  Dept.. 

Iowa  State  University 
Amrs,  Iowa  50010 

Aorcspacw  Engiiiearing 
Dept. 


I 

I 

i 

! 

} 

I 

) 

I 


I 


! 


1 


I 

( 

I 

: 

I 

9 


i 


i 


D1STM8UT20W  (COW) 

Copies 


The  Johns  Hopkins  University 
Baltimore,  Maryland  21218 
Prof.  s.  Corroin 

University  of  Kentucky 
Kenner-Cron  Aero.  Lab. 

Lexington,  Kentucky  40506 
C.  r.  Knapp 

Department  of  Aero. 

Bngineering,  ME  106 
U)uisianna  State  University 
Baton  Rouge 
Louisiana  70303 
Or.  P.  H.  Miller 

University  of  Maryland 
College  Park 
Maryland  20740 

Prof.  A.  Wiley  Shervood 
Oepartment  of  Aerospace 
Engineering 

Prof.  Charles  A.  Shreeve 
Department  of 
Mechanical  Engineering 
Or.  S.  1.  Pai,  Institute 
for  Fluid  Dynamics  and 
Applied  Mathematics 
Dr.  Redfleld  M.  Allen 
Departnent  of 
Mechanical  Engineering 
Dr.  W.  L.  Melnik 
Oepartment  of 
Aerospace  Engineering 
Dr.  John  D.  Anderson,  Jr, 
Oepartment  of 
Aercspace  Engineering 

Michigan  State  University 
Library 
East  Lansing 
Michigan  48823 
Documents  Department 

Ksssachusctts  institute  of 
Technology 
Casdjridge 

Massachusetts  02139 
Mr.  J.  R.  Hartuccelli 
Rm.  33-211 
Prof.  M.  Finston 
Prof.  J.  Baron,  Dept, 
of  Aero,  and  Astro. 

»a.  37-461 
Prof.  A.  H.  Chapiro 
Itsad,  Mach.  Cngr.  Dept. 

Aero.  Engineering  Library 
Prof.  Ronald  Probestein 
Or.  E.  E.  Cc/ert 

Aerorhysics  Laboratory 

Daiversity  of  Michigan 
Ann  Arbor,  Michigan  48104 

Dj . M.  sickel.  Dept  ot  Aero  Engr 
Engineering  Library 
Aerospace  Engineering  Lib. 

Mr.  C.  CousincaUr  Engin-Trans  Lib. 


Serials  and  Oocumrnte 
Section 

Cenercl  Library 
University  of  Michigan 
Ann  Arbor,  Michigan  48104 

Kiaaissjppi  Stetu 
vinivorsity 

Oepoitment  of  Aercphyaics 
and  .lororp.ice  Engineefing 
P.O.  Ora»vr  A 

ftate  College,  Misslnsippi  39762 
Mr,  n.  cijctt 


U.S.  Naval  Academy 
Annapolis,  Ma.-yland  21402 
Engineering  Department 
Aerospace  Divis'on 

Library,  Code  2124 
U.  S.  Naval  Postgraduate 
School 

Monterey,  California  93940 
Technical  Reports  Section 

New  York  University 
University  Heights 
New  York,  Hew  York  10453 
Dr.  Antonio  Fcrri 
Director  of  Guggenheim 
Aerospace  Laboratories 
Prof.  V,  Zakkiy 
Engineering  and  Science 
Library 

Morth  Carolina  State  College 
Raleigh 

North  Caiolina  27607 

Or.  F.  K.  OeJarnette,  Dept 
Mech.  and  Aero. 
Engineering 

Dr.  H.  A.  Hassan,  Dept,  of 
Mach,  and  Aero.  Engr. 

D.  H.  Hill  Library 
Morth  Carolina  State 
University 
P.O.  Box  5007 
Raleigh 

North  Csrolins  27607 

University  of  North  Csroli"a 
Chapel  Hill 
Morth  Carolina  27514 
Oepartment  of  Aero. 
Engineering 

Library,  Documents  Section 
AFROTC  net  390 

Northwestern  University 
Technological  Institute 
Bvanr.ton,  Illinois  60201 
Department  of  Hechonica: 
Englnesring 
Library 


Virginia  Polytechnical  Instltuta 
Blacksburg,  Vc.  24061 
Prof.  G,  Inger 


Department  of  Aero -Astro 
Engineering 
Ohio  Stats  '.nivereity 
2036  Neil  Avenuj 
Colu.mbu3,  Ohio  4J210 
Engineer  <r.g  Library 
Prof.  J.  D.  Lee 
Prof.  G.  L.  Voi-  Etcher. 

Stete  University 
Libraries 
Documents  Divis'on 
1858  Noli  Avo'ue 
Columbus,  Ohi.  43210 

The  Pcimsylvania  Stote 
University 
University  Part 
Penney  Ivan.’ a 18602 
Dept-  of  Aaro  Engr. 

Hrmmond  Bldy. 

Librsr/j  Documents 
Bastion 

Bsvier  Eitginenrin.j  Library 
135  Oc-redum  Hall 
University  cf  PSttrb.'rgh 
Pittsburgh 
Pennsylvsnia  15262 


2 


DibrnmurioN  (cont) 


COPIES 


COPIES 


r 


I 


> 

I 

c 

i 

r 

? 

5 

I 


i 


i 


' '.k'elo.'i  ^ '\f  'c  v'y 

.... ( SciencQ  Ccpt. 
D-2.  ! l.iKii.j.  ‘...i.ulranqle 
Prii  .cton 
New  .Jersey  08S40 
Prof.  S.  Bogdonoff 
Dr.  I.  E.  Vas 

Purdue  University 
School  of  Aeronautical  -and 
I'jiqineering  Sciences 
Lafayette,  Indiana  47907 
Li brary 

Dr.  B.  Rvese,  lic.id.  Dept 
of  Aero.  * Astro. 

Rensselaer  Polytechnic 
Institute 

Troy,  New  York  12181 
Dept,  of  Aeronautical 
Fngincci'inq  and 
Astronautics 

Dcpartiaent  of  Mechanical 
Industrial  and  Aerospace 
Engineering 
Rutgers  - The  State 
Dniversity 

New  Brunswick,  ;i.  J,  08903 
Dr.  R.  U.  Page 
Dr,  C.  p.  Cnen 

Stanford  University 
Stanford 

California  9430S 
Librarian,  Dept,  of 
Aeronautics  and 
Astronautics 

Stevens  Institute  of 
Technology 

Boboken,  New  Jersey  07030 
tiechanical  Engineering 
Department 
Library 

The  University  of  Texas 
at  Austin 

Applied  Research  Laboratories 
P.  0.  Box  8029 
Austin,  Texas  78712 
Director 

Engr  S.B.114B/Dr.  Friedrich 

University  of  Toledo 
2S01  k.  Bancroft 
Toledo,  Ohio  43606 
Dept,  of  Aero 
Engineering 
Dep*.  of  Mcch 
Engineering 

University  of  Virginia 
Ecboc/1  of  Ergineeting 
and  Applied  Science 
Charlottesville 
Virginia  22901 
Dr.  I.  D.  Jacobson 
Dr,  G.  Matthews 
Dr.  R,  H.  Zapat.i 
University  of  Washington 
Seartls 

Waahlngton  98106 
Engineering  Library 
Dept,  of  Aeronautics  and 
Astron.actics 

prof.  R.  E.  .street.  Dept, 
of  Aoro.  and  Astro. 

Prof.  A.  .'tsrtzborg.  Aero, 
snd  Astro. , Cuggnhcin 
llsll 

Nest  Vlr.ijnfa  University 
Mergantowr. 

Nest  Virainia  2C606 
Library 


FedernJ  Reports  Center 
University  of  Wisconsin 
Hcch.i:iicnl  Engineering 
Building 

Madison,  Wisconsii.  53706 
S,  Reilly 


Prototype  Development  Associates 
1740  Garry  Avenue 
Suite  201 

Santa  Ana,  C.\  9270S 

Dr.  J.  Dunn 
Dr.  P.  Crenshaw 


Los  Alamos  Scientific 
I.aborotory 
P.O.  Box  1663 
Los  Alamos 
New  Mexico  87544 
Report  Library 

University  of  Maryland 
Baltimore  county  (UMBO 
5401  Nilkens  Avenue 
Beltinorc,  Maryland  21228 
Dr.  R.  C.  Roberts 

Mathematics  Departsient 

Systems  Research  Laboratories,  inc. 
3800  Indian  Ripple  Road 
Dayton,  Ohio  4S440 
Dr.  X.  Ball 
Dr.  C.  Ingram 

Institute  for  Defense 
Analyses 

400  Army-Navy  Drive 
Arlington,  Virginia  22202 
Classified  Library 

Xanan  Sciences  Corporation 
P.O.  Box  7463 
Colorado  Springs 
Colorado  80933 
Library 

Xaman  Science  Corporation 
Avidyne  Division 
S3  Second  Avenue 
Burlington 

Ksssachusetts  01803 
Dr.  J.  R.  Ruetenik 

Rockwell  International 
B-1  Division 

Technical  Information  Center 
(BA08) 

International  Airport 
Los  Ann'-les,  Ca.  90009 

Rockwell  International 
Corporation 

Technical  Information  Center 
4300  E.  Filth  Avenue 
Columbus,  Ohio  43216 

M.  1.  T.  Lincoln  Laboratory 

P.O.  Box  73 

Lexington 

Massachusetts  02173 
Library  A-082 

The  RAND  Corporation 
1700  Main  Street 
Santa  Monica 
California  99406 
Library  - D 

Aorojet  Eloctrosystc-es  Co. 

1100  1*.  Eij]  lyv.ilu  Avu. 

Axusu,  Ca.  91702 
Engineering  Library 


1 


m 


c.  . 


DJSTnillUTION  (CONT) 


i 


i 


1 


Copiea 


Copl«a 


The  Boeing  Company 
P.O.  Box  3999 

Seattle,  Waaliington  98124 
87-C7 


Unitud  Aircraft 
ftCBoarch  Laboratoriea 
East  Hartford 
Connecticut  0G108 
Dr.  William  M.  Foley 

U.iited  Aircraft  Corporation 
400  Kain  Street 
w«st  Hartford 
Connecticut  06108 
Eihrary 

Hughes  Aircraft  Company 
Ccntincla  at  Toale 
Culver  City,  Ca.  90230 
COT.pany  Tech.  Doc.  Center 
S/Ell,  B.  W.  Campbell 

Lockheed  Missiles  t Space  Co. 
Continental  Bldg.,  Suite  445 
El  Segundc,  CA  90245 
T.  R.  Fortune 
F.  £.  Huggin 


Lockheed  Missiles  and  Space 
Company 
P.O,  Box  504 
Sunnyvale 
California  9408C 

Hr.  G.  M.  Laden,  Dept. 

81-25,  Bldg.  154 
Kr.  M:]rl  Culp 

Lockheed  Missiles  and  Space 
Company 

3251  Hanover  Street 
Palo  Alto,  California  94304 
Technical  Information 
Center 

Lockliecd-Cal  i fornia 
Company 

Burbank,  California  91503 
Central  Library,  Dept. 
•4-40,  Eldg.  170 
PLT.  8-1 

Vice  President  and  Chief 
Scientist 
Dept.  03-10 
Lockheed  Aircraft 
Corporation 
P.O.  :*ox  551 

Burbank,  California  91503 

Hartlr  Marietta  Corporation 
P.O.  Box  988 
Baltlmoro 
Maryland  21103 

Science-Technology  Library 
(Mail  Ho.  398) 

Martin  Company 
3211  Trade  Winds  Trail 
Orlando,  F.lorida  32805 
Mr.  1!.  J.  Dicbolt 

General  Dynaiaies 
P.O.  Box  748 
Fort  Worth,  Texas  76101 
Research  Library  224C 
Ceoree  Xaler,  Mail  Zone 
2810 


Calsp.in  Corporation 
4455  Genesee  Street 
Buffalo,  New  York  14221 
Library 

Air  Univeraity  Library 
(SE!  63-578 

Maxwell  Air  Force  Bate 
Alabama  36112 

McDonnell  Company 
P.O.  Box  516 

St.  Louis,  Missouri  63166 
R.  D.  Dctrich,  Dept.  209 
Bldg.  33 

W.  Brian  Broo)cs 

McDonnell  Douglas  Astronautics  Co.  - West 
5301  Bolca  Avenue 

Huntington  Beach,  California  92647 
A3-339  Library 
3.  S.  Murphy,  A3-B33 
M.  Michael  Briggs 


Peirchlld  Hiller 
Republic  Aviation  Diviaion 
Panaingdale 
New  York  11735 

Engineering  Library 

General  Applied  Science 
Laboratoriea,  Inc. 

Marrick  and  Stewart 
Avenues 

Nestbury,  Long  island 
Maw  York  11590 
Or.  P.  Lane 

L.  M.  Nucci 

General  Electric  Company 
Research  and  Develcp.ment 
Lab.  (Comb.  Bldg.) 

Schenectady 
Maw  York  12301 

Or.  H.  T.  Nagamatsu 

The  Whitney  Library 
General  Electric  Reaearch 
and  Development  Center 
..  Knolla,  K-1 
P.O.  Box  8 
Schenectady 
Mew  Ycrk  12301 

M.  P.  Orr,  Manager 

General  Electric  Company 
Miatile  and  Space  Division 
P.O.  Box  8555 
Philadelphia 
Pennsylvania  19101 
MSD  Library 

Larry  Chaien,  Mgr. 

Or.  J.  D.  Stowart,  Mgr. 
Research  and  Engineers 7g 

General  Electric  Company 
AEG  Technical  Information 
Center,  M-32 
Cincinnati,  Ohio  45215 

General  Electric  Company 
Re-Entry  i Environmental  Systema 
••ivlsicn 

Chestnut  Street 
Philadelphia,  Punn.  19101 
Or.  S.  M.  Seals 
Dr.  II.  Lew 
Hr.  J.  W.  Pauat 
A,  Hartelluscl 
W.  Daskin 
J.  D.  CrcsBwcll 
3.  PottUN 
L.  A.  .M<,rbhall 
3.  C.iRiBnto 
R.  l(o!>l.g: 

C.  IMrrie 
F.  Cc-orge 


! 

i 

5 


4 


DisntiBUTXON  (ccart) 


I 


AVCO-Ev«rett  Rrgecrch 
Lftboratcry 

'/38S  Revera  B««eh  Parkway 
Everett 

Mciaechueettp  0314$ 

Library 

Or.  Ceor^e  Sutton 

LTV  Aeroapaee  Corporation 
Vougbt  Eeronautlca  Olviklon 
P.O.  Box  S907 
Dallsa,  Texas  75332 
Onlt  2-51131  (Library) 

LTV  Aerospace  Corporation 
Mlsoiles  and  Space  Division 
P.O.  Box  32S7 
Dallas,  Texas  75223 
KSD-T-Library 

Northrop  Norair 
3501  West  Broadway 
j Hawthorne 

California  $0250 
Tech.  Info.  3360-32 

CoverniMsnt  Documents 
The  Pcundren  Library 
Nice  Institute 
P.O.  Box  1«92 
Booston,  Texas  77001 

CrusoMn  Aircraft 
Engineering  Corporation 
Bethpaga,  Long  Island 
Mew  York  11714 
Mr.  B.  A.  Sctieuing 
Mr.  B.  B.  Hopkins 
Mr.  B.  X.  Xasd 

Mar<iuardt  Aircraft 
Corporation 
1C55S  Saticoy  Street 
Van  Kuys,  California  91409 
Library 

AXDE  Associates 
P.O.  Box  286 
510  Mlnters  Avenue 
Poranus,  New  Jersey  07653 
Librarian 

Aerophyslcs  Coaipany 
3500  Connecticut  Ave.,  N.M. 
Mashlngton,  D.C.  20003 
Mr.  C.  0.  Boehler 

Aeronautical  Research 
Associates  of  Princeton 
50  Washington  Road 
Princaton 

Mew  Jersey  01540 

Dr.  C.  duP.  Donaldson 

General  Research  Corporation 
5313  Hollister  Avenue 
P.O.  Box  3587 
Santa  Barbara 
California  93105 
Technical  Infonwtlon 
Office 

Sandia  Laboratories 
, Nall  Service  Section 

• Albixjuorque,  3,  h.  87115 

Mr.  K.  Coin,  Div.  5262 
Mr.  M.  H.  Curry,  Olv.  1331 
Hr,  A.  M.  Tornsby,  3141 
Or.  C.  stone 
Div,  3141 

Hercules  Incorporated 
Allegany  BalHstlcs 
Laboratory 
P.O.  Box  210 
Cuaborlsnd 
Maryland  31503 
Library 


Copies  Copies 

General  Electric  Coapany 
P.O.  Box  2500 
Daytona  Beach 
Florida  32015 

Dave  Hovis,  Ka.  4109 

2 

nil  Systems  Group 
1 Space  Park 
Redondo  Beach 
California  90278 

Technical  Libr//Doc  Acquisitions 
B.  Pearce,  Aero  Dept. 

P.  D.  Deffanbaugh 
Stanford  Research  Institute 

333  Ravsnswood  Avenue  ‘ 

Menlo  Park  ■ 

California  94025  I 

Dr.  C.  Abrahamson  i 

Hughes  Aircraft  Company 
P.O.  Box  3310 
Fullerton 
California  92634 

Technical  Library,  600-C232 

Mastinghouse  Electric 
Corporation 

Astronuclear  Laboratory 
P.O.  Box  10864 
Pittsburgh 
Pennsylvania  15236 
Library 

Oniverslty  of  Tennesse 
Bpaee  Institute 
Tullahona 
Tennessee  37388 
.Prof.  J.  M.  Hu 

COHVAIR  Division  of 
General  Dynsnics 
Library  and  Infonsation 
Services 
P.O.  Box  12009 
San  Diego 
California  92112 

CCMVAIR  Division  of 
General  Dynsaics 
Post  Office  Box  0986 
San  Diego,  California  92138 
Dr.  J.  Raat 
Mail  lone  640-02 
Research  Library 

AVCO  Missiles  Systens 

' Division  - i 

301  Lowell  Street 
W1  Islington 
Massachusetts  01887 
E.  E.  H.  Schursana 
J.  Otis 

Chrysler  Corporation 
Space  Division 
P.O.  Box  2920C 
New  Orleans,  La.  70189 
H.  D.  Kenp,  Dept.  2910 
E.  A.  Rawls,  Dept.  2920 

General  Dynamics 
PosK>na  Divisicn 
P.O.  Box  3507 
Pomona,  Ca.  91766 
Tech.  Doc.  Center,  Mall  xone  6-20 

General  Electric  Cos«>aay 
3196  Chasnut  Street 
Philadelphia,  Pa.  19101 
H.  Danskin 
Larry  Chasen 

Dlte  He  XfCV 

Phllco-Ford  Corporation 
Aaroneutronic  Division 
Newport  Beach 
California  92660 
Dr.  A.  Demotriadss 


1 

f 


J 


1 


DISTRIBUTION  (COOT) 


i 


■ i 
' i 

t 


Copiaa 

Raytheon  Company 
Missile  Systems  Division 
Hartwell  Road 
Bedford,  Ma.  01730 
D.  P.  Forsmo 

TRW  Systens  Group 
Space  Park  Drive 
Houston,  Texas  770S8 
M.  H.  Sweeney,  Jr. 

Marine  Bioscience  I.aboratory 
513  Sydnor  Street 
Rid9ecrc8t,  Ca.  93555 
Dr.  A.  C.  Charters 

University  of  California  - 
Los  Angeles 

Dept  of  Mechanics  t Structures 
Los  Angeles,  Ca.  90024 
Prof.  J.  D.  Cola 

University  of  Hyomins 
University  Station 
P.  O.  Box  3295 
Laramie,  Wyoming  82070 
Head,  Dept.  Mech.  Cng. 

Applied  Mechanics  Review 
Southwest  Resesrch  Institute 
8500  Culebra  Road 
San  Antonio,  Texas  78228 

American  Institute  of  Aeronautics 
and  Astronautics 
1290  Sixth  Avenue 
Hew  York,  Hew  York  10019 
J.  Newbauer 

Technical  Information  Service 
American  Institute  of  Aeronautics 
and  Astronautics 
750  Third  Avenue 
Hew  York,  Hew  York  10017 
Miss  p.  Marshall 
Faculty  of  Aeronautical 
Systems 

University  of  West  Florida 
Pensacola,  Florida  32504 
Or.  R.  Fledderman 

Space  Research  Corporation 
Chittenden  Bank  Building 
North  Troy,  Venaont  05859 
Library 
J.  A.  Finkal 

««e  Aerospace  Corporation 
P.  O.  Box  92957 
Xx>s  Angeles,  Californio  90009 
J.  M.  Lyons,  Bldg.  82 

Chrysler  Corp. , Defense  Division 
Detroit,  Michigan  48231 
Or . R.  Lusardi 


AERO 

3020  BueklnghsM  Drive 
South  Bend,  Indians  4(814 
Dr.  J.  Nicolaidan 

Acurax  Corp.  Aerothans 
485  Clyde  .Aveuua 
Mt.  View,  CA  94042 
L.  Cooper 

Sandla  Corporation 
Livermore,  CA  94550 
J.  K.  Xryvoruka 


Copies 

Near,  Xnc. 

510  Clyde  Avenue 
Mountain  View,  CA  94043 


^09^AIR  Division  of  General  Dynasties 
P.O.  Box  80847 
San  Diego,  California  92138 
Dr.  Z.  S.  Levinaky 
Kail  Zona  (S7-1 


.■■'fern  »■>  .'ll  ji  I nwKseeiii  .jim- 


4 


